Topological insulators and metal-insulator transition in the pyrochlore iridates 
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The possible existence of topological insulators in cubic pyrochlore iridates A2lr207 (A = Y or rare-earth 
elements) is investigated by taking into account the strong spin-orbit coupling and trigonal crystal field effect. 
It is found that the trigonal crystal field effect, which is always present in real systems, may destabilize the 
topological insulator proposed for the ideal cubic crystal field, leading to a metallic ground state. Thus the 
trigonal crystal field is an important control parameter for the metal-insulator changeover. We propose that this 
could be one of the reasons why distinct low temperature ground states may arise for the pyrochlore iridates 
with different A-site ions. On the other hand, examining the electron-lattice coupling, we find that softening 
of the q=0 modes corresponding to trigonal or tetragonal distortions of the Ir pyrochlore lattice leads to the 
resurrection of the strong topological insulator. Thus, in principle, a finite temperature transition to a low- 
temperature topological insulator can occur via structural changes. We also suggest that the application of the 
external pressure along [111] or its equivalent directions would be the most efficient way of generating strong 
topological insulators in pyrochlore iridates. 
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I. INTRODUCTION 

Topological band insulators arise from nontrivial Berry 
phase of electron wave functions and possess gapless bound- 
ary states as a consequence of topological properties of the 
bulk electron energy bands. '"'^ Topological invariants of elec- 
tron wave functions in such topological insulators can be used 
to describe/identify two dimensional quantum spin Hall in- 
sulators and three dimensional strong topological insulators. 
The presence of the gapless boundary states and the associ- 
ated topological properties have recently been confirmed by 
a series of remarkable experiments on HgCdTe, Bii-j^Sb^;, 
Bi2Se3, Bi2Te3, and other materialsJ^— It has been known 
that the strong spin-orbit coupling in these systems provides 
an essential ingredient for the nontrivial Berry phase of elec- 
tron wave functions. One important future direction is to un- 
derstand the effect of interactions on the topological insula- 
tors. Growing interests on combined effects of electron corre- 
lation and topological properties have naturally lead to recent 
fascinating studies on the realization of topological phases in 
transition metal oxides with d electrons,''''^" instead of more 
conventional s or p orbital systems where electron coiTelation 
is less important. 

Transition metal oxides with 5d electrons are characterized 
by the strong spin-orbit coupling due to the large atomic num- 
ber of 5rf transition metal elements. As a result, the spin-orbit 
coupling competes with the kinetic and interaction energies, 
leading to substantial correlation effects despite the relatively 
extended nature of the 5d orbitals. This feature, for example, 
is confirmed by recent experiments on Sr2lr04^^^ where the 
spin-orbit coupling plays an essential role in the formation of 
the Mott-insulator ground state. In the case of the so-called 
hyperkagome material Na4lr308, a spin liquid ground state 
is proposed, manifesting strong coiTelation effects in 5d elec- 
tron sy stems. ^^"^^^ Moreover iridium oxides are considered as 
promising candidate systems where we can study the inter- 
play between electron correlation and strong spin-orbit cou- 
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FIG. 1 : (a) Structure of the pyrochlore lattice. An Ir atom sits at each 
vertex of a tetrahedron. Four sites within a tetrahedral unit cell are 
numbered from 1 to 4. (b) The coordinate system used in this paper. 



pling. Perhaps most interestingly, the iridium oxides may also 
be ideal materials for the occuiTence of topological insula- 
tors. It is suggested that the honeycomb lattice of Ir ions in 
Na2lr03 with the complex hopping amplitudes arising from 
the strong spin-orbit coupling may lead to a quantum spin 
Hall insuIatorJii Three dimensional pyrochlore lattice of Ir is 
also suggested as a candidate system for a strong topological 
insulator.^" The possibility of the topological Mott insulators 
where the spinons, not the electrons, possess topological band 
structures is discussed as a result of correlation effectSi^Si^ 

In this work, we study the possibility and stability of strong 
topological insulators in pyrochlore iridates, A2lr207 (A = 
Y or rare-earth elements) by taking into account both of the 
strong spin-orbit coupling and trigonal crystal field effect in- 
duced by a local distortion of IrOg octahedra. Due to the ex- 
tended nature of 5d electron wave functions, the local crys- 
tal field strength in 5d systems also becomes a large energy 
scale competing with electron interaction and spin-orbit cou- 
pling. The strength of the trigonal crystal field effect may 
be different for distinct choices of the A-site ion in A2lr207. 
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Given that there exist a variety of choices for the A-site, the 
pyrochlore mdates offer the opportunity to control the relative 
strength of the crystal field effect and the spin-orbit coupling. 
In addition, the pyrochlore iridates show both metallic and 
insulating ground states at low temperatures as well as finite 
temperature metal-insulator transitions depending on the A- 
site ions.^^"^^ Thus it would be very interesting to understand 
the origin of distinct phases and phase transitions. 

Here we derive an effective Hamiltonian for the Ir ions in 
the pyrochlore iridates by carefully incorporating the trigonal 
crystal field effect and using the appropriate spin-orbital ba- 
sis. Previously an effective Hamiltonian was derived by Pesin 
and Balents^" for the Ir ions in a perfect octahedral environ- 
ment with the cubic Oh symmetry and it was found that the 
Ir pyrochlore system supports a strong topological insulator 
when the strength of spin-orbit coupling is sufficiently large. 
In contrast we find that the trigonal crystal field effect, which 
is always present in real pyrochlore systems, is quite signifi- 
cant and it can lead to a metallic ground state. 

On the other hand, we find that the distortions of the Ir py- 
rochlore lattice induced by softening of certain q=0 phonon 
modes give rise to the resurrection of strong topological in- 
sulators. This suggests that finite temperature metal-insuator 
transition can occur through structural distortions such that the 
low temperature ground state is a strong topological insulator. 
It is interesting to notice that recent experiments on Sm2lr207 
and Eu2lr207 show structural changes at the metal-insulator 
transition even though the low temperature ground state seems 
to be magnetic and show spin-glass-like behavior ^^•^'^ We 
also suggest from the studies of the electron-lattice coupling 
that applying pressure along [111] or its equivalent directions 
would be the most efficient way of generating strong topolog- 
ical insulators in pyrochlore iridates. 

The rest of the paper is organized as follows. In Sec.HIlwe 
derive the effective hopping Hamiltonian including the trigo- 
nal crystal field effect. We first consider the influence of local 
trigonal distortion of the oxygen octahedra on Ir <2g electrons. 
Then we explain the derivation of the corresponding lattice 
Hamiltonian. In Sec.|lIIl the evolution of the electronic struc- 
ture induced by trigonal crystal field in the presence of the 
strong spin-orbit coupling is described in detail. Here we ex- 
plain how a metallic state may arise when the trigonal crystal 
field effect is present. In Sec. IIVI we discuss the effect of 
electron-lattice coupling on the electronic structure of the Ir 
pyrochlore system. We show how the electron-lattice coupling 
can lead to strong topological insulators. Finally, we conclude 
in SeclV] The detailed expressions for various matrices de- 
scribing the effective Hamiltonian are given in the Appendix. 



II. EFFECTIVE HOPPING HAMILTONIAN WITH 
TRIGONAL CRYSTAL FIELD 

In this section we investigate the effect of local trigonal 
crystal field on the electronic structure of Ir 5d electrons. In 
Sec. Ill Al we explain the degeneracy lifting of t2g orbitals un- 
der trigonal crystal fields. The effective hopping Hamiltonian 
including the trigonal crystal field is derived in Sec. Ill Bl De- 




(b) (c) 




FIG. 2: Trigonal crystal field around an Ir ion. (a) Local geome- 
try of an IrOe octahedron. A trigonal distortion is induced by com- 
pression or elongation of the surrounding oxygen octahedron along 
a C-j. symmetry axis for 3-fold rotation, (b) Splitting of degenerate 
t2g States under the trigonal crystal field, (c) Local C3 axes for four 
Ir ions constituting a tetrahedron; Each C3 axis points to the body 
center of the tetrahedron. 



tailed procedures of the derivation can be found there, as well. 



A. Trigonal crystal field 

Each Ir^"*" ion is coordinated by six oxygen anions which 
are at equal distance from the central Ir**^ cation. The actual 
coordinates of the surrounding oxygen anions have a free po- 
sitional parameter, so-called the oxygen x parameter, which 
depends on material properties. In general, the six oxygen 
ions around a central Ir^+ form a distorted octahedron where 
the amount of distortion can be parameterized by the oxygen 
X parameter. For x — Xc = 5/16, each Ir*^ ion is under a 
perfect local cubic crystal field. The deviation of x from the 
ideal value of Xc generates a trigonal crystal field. The trig- 
onal distortion of an oxygen octahedron is accompanied by 
compression {x > Xc) or elongation {x < Xc) of the oxygen 
octahedron along one of the C3 symmetry axis for 3-fold ro- 
tation. A local geometry of an IrOe octahedron and a C3 sym- 
metry axis for the trigonal distortion are described in Fig.|2ja). 

Under the trigonal distortion, the local symmetry group 
around an Ir site reduces from the cubic Oh to D^d- Due 
to the trigonal crystal field, degenerate t2g orbitals split into 
a singlet aig state and a doublet e'g state belonging to Aig 
and Eg irreducible representations of the D^d point group, re- 
spectively. (See Fig.|2|b).) A C3 symmetry axis corresponds 
to one of the [111] or its equivalent axes. Taking the local 
[nin2n3] directions as the C3 symmetry axes (n^ = ±1), the 
effect of trigonal crystal field on t2g states can be described 
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by the following Hamiltonian, 




nin2 riin^ 
711^2 n2n3 
nins n2n3 



Here we have chosen a local basis as 5*^ = {d}y^, d\^, d^y). 
The singlet aig state with the energy Ea — — 2Atji/3 is given 
by 
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|ai) = al|0) = -j=[ni\dy^) + n2|4a;) + "-sMa^y)]- (1) 

The eigenvectors for the e'g doublet states with the energy 
Ee — Atri/3 can be written as 
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|e+) = e^lO) = -j=[niLo\dy:,) + n2i^ + "-sMxy)], 

|e_) = el|0) = ■^[niuj'^\dy^) + n2Uj\dzx) + nz\dxy)], 



(2) 
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FIG. 3: Local Ir-O-Ir configuration. Local octahedral axes for Ir 
atoms on the left-hand side (Iri) and the right-hand side (Irn) are 
not parallel to the global cubic axis adopted for the central oxygen 
atom. 



around an Ir site without a trigonal crystal field, ta = is sat- 
isfied. In other words, the effect of trigonal crystal fields on 
the electronic structure can be captured by varying the relative 
magnitude of ia/^e- In the following section, we explain the 
procedures deriving the above effective Hamiltonian in Eq.(|3]|. 



where lo = e*~. 

In Fig. EJc), we show the local C3 axes of the four Ir sites 
constituting a tetrahedral unit cell. Each C3 axis points to the 
body center of the tetrahedron. Considering the trigonal dis- 
tortions of IrOg octahedra and their relative orientations, we 
construct the effective hopping Hamiltonian for Ir d electrons 
on the pyrochlore lattice. The resulting Hamiltonian is given 
by 



(3) 



B. Construction of the effective hopping Hamiltonian 

We construct the effective Hamiltonian describing hop- 
ping processes between neighboring Ir t2g orbitals through 
intermediating oxygen p orbitals in the following way. We 
first consider a small cluster composed of neighboring Ir-O- 
Ir atoms as shown in Fig. [3] Allowing Ir-O hopping pro- 
cesses, the three-atom cluster can be described by the follow- 
ing Hamiltonian, 



(ij) n,n' a, a' 

where 

Here i and j are indices for unit cells while n and n' indi- 
cate Ir sites within a unit cell (n — 1,2,3,4). We use the 
index a (a — 1, ...6) to describe local spin-orbit eigenstates. 
Within the t2g manifold, d electrons behave as if their ef- and 
fective orbital angular momentum is one with an additional 



minus sign.^" Therefore the effective total angular momen- 
tum of Ir d electrons, jeff can be either 1/2 or 3/2. We use 
a — 1,2 to indicate the spin-orbit doublet with jeff — 1/2 
and a — 3,4, 5,6 to denote the spin-orbit quadruplet with 
jeff = 3/2. The on-site energy = Xso for the spin-orbit 
doublet and = —Xso/"^ for the spin-orbit quadruplet with 
the spin-orbit coupling strength given by \so- Tna,n'a' is the 
hopping amplitudes between nearest neighbor Ir sites and /i is 
the chemical potential. 

In contrast to the tight-binding Hamiltonian derived by 
Pesin and Balents in Ref20 where a single energy scale 
t describes the hopping processes, the hopping amplitude 
Tna,n'a' has two independent parameters ta and in this 
case. These two hopping parameters come from the hopping 
processes between aig and states under the trigonal crystal 
field, respectively. Under the perfect local cubic environment 



H^H^°'>+H^^\ (4) 



= tjIVlP + P^VlvL + 7]^iiVRp + p^Vl'nB. (6) 

Here H'-^'> indicates on-site potentials at Ir and O sites while 
H'-^'f describes Ir-O hybridizations. The on-site atomic po- 
tentials of d and p orbitals give rise to the diagonal matri- 
ces, H^°'> = dmg[Ea,Ee,E,] and H^°'> = diag[Ep, Ep, Ep]. 
Here Ea (Ee) is the atomic energy for an aig singlet (e'g dou- 
blet) state composed of Ir d orbitals under the trigonal crys- 
tal field and Ep is the atomic energy of degenerate oxygen 
p orbitals. Vn (n= L, R) represents overlap integrals be- 
tween Ir d orbitals and oxygen p orbitals. Taking into ac- 
count the local trigonal crystal field of each Ir atom, we em- 
ploy a basis, that is diagonal under the local trigonal crys- 
tal field, ry^ = {a\ ,^,e\^ el_ „) to represent Ir t2g orbitals. 
The subscript n indicates the locations of two Ir atoms on 
the left-hand side (n = L) or the right-hand side (n = R). 
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P' — (Px T Py ' pI) represents oxygen p orbitals connecting 
neighboring Ir atoms. It is important to notice that the local 
octahedral axes of Ir atoms are not parallel to the global cubic 
axes as described in Fig. [3] Each component of the oxygen p 
orbitals is defined with respect to the global cubic axis, while 
d orbital basis is defined with respect to the local octahedral 
coordinate of each Ir atom. 

The effective hopping Hamiltonian between two neighbor- 
ing Ir atoms is obtained from the second order perturbation 
theory treating H'^-^^ as a perturbation. The resulting Hamilto- 
nian is given by 



(0) 



H 



(0) 



(0) 



H 



(0) 



|-7/_R + h.c. 

(7) 



The information about the relative orientation between two 
local octahedral axes of neighboring Ir atoms is contained in 
Vl and Vr. To construct Vl and Vr, we follow the approach 
used by Pesin and Balents in Ref. 20. We first construct ma- 
trices describing hopping processes between Ir t2g and O p 
orbitals. It is convenient to define components of oxygen p 
and Ir t2g orbitals with respect to the same cubic axis. In 
this case the Ir-O hopping processes are strongly constrained 
by the symmetries of d and p orbitals. The hopping matrices 
connecting {p^, Py, Pz) orbitals to {dyz, dzx, dxy) orbitals are 
given by^ 




pdTT 





dxy) to the trigonal basis (ai, e+, eJ) in the following way, 
From Eq. ([T]i and Eq. (|2|i, is obtained as. 



n d'^ 



^ \ nicj^ n2U! 



Now we can construct Vl and Vr explicitly, which are given 
by Vl = 17*rp(i)i?(^)t and Vr = n*Tp^R^R'-"^^l Here the 
subscript j'(L) (p{R)) of r indicates the location of the inter- 
mediating oxygen atom with respect to the local octahedral 
axis of the Ir ion on the left-hand (right-hand) side. 

The effective hopping Hamiltonian in Eq.© can be written 
compactly as. 



Hlm = lvl{ThopVLVl + 14l4^hop}??i?^ + h., 



(8) 



where Thop=diag[ta,te,te] with ta 



:JiEa - Ep) and 



''pd 



-^/{Ee — Ep). In the above, we have replaced Vr 
{Vl) by VRjapdn {VL/apdn) to make ta and tg have conven- 
tional forms for hopping amplitudes. Note that ta = t^ in the 
absence of the trigonal crystal field effects because Ea = E^ 
is satisfied in that case. The effect of local trigonal field split- 
ting can be considered by changing the relative magnitude of 

Up to now we have neglected spin degrees of freedom. 
Since the spin-orbit coupling is the largest energy scale of the 
problem, it is convenient to use the local spin-orbit eigenstates 
as a basis for the representation of the Hamiltonian. To project 
the Hamiltonian onto the local spin-orbit basis, we define ro- 
tation matrices Z?*^*' acting on the spin space, which is noth- 
ing but a spinor representation of 0(3) rotations i?*^*' defined 
before. The detailed expressions of matrices i?*^*' and D^^^ 
are given in the Appendix. Then the spin dependent hopping 
Hamiltonian is given by 



Hl,r^ ^ v[^ia'^i<T,i''7'VR,i'<7' + h.c, (9) 



lcr,l' 



where 



Here the subscript and superscript of the matrix r indicate 
the location of an O atom with respect to the Ir atom. For ex- 
ample, r+ describes Ir-O hopping processes when an O atom 
sits on the positive x axis with the Ir atom sitting at the ori- 
gin. apd-K is the hopping integral between neighboring p and 
d orbitals. 

Since the local octahedral axes of Ir atoms are not paral- 
lel, we introduce an additional matrix i?*^*^ for each ith Ir 
atom, which rotates the global cubic axis to a local octahe- 
dral axis. Under the action of i?^*-', the pn (n = x, y, z) or- 
bitals defined with respect to the global cubic axis transform 



as p. 



Rmnpln- Here Pi n denotes a p orbital defined with 
respect to the local octahedral axis of the ith Ir atom. In ad- 
dition, we introduce the matrix 51 that transforms {dyz, dzx, 



T, 



1(7.1' I 



-,[n,pVLvl + VLV^n,^]u,[{D^'^^)^D^% 



(10) 



Here I is the index for d-orbitals under the local trigonal crys- 
tal field (l = a, e+, e_) and a indicates the spin projection 
with respect to a local quantization axis. 

We can obtain the final expression for the effective hop- 
ping Hamiltonian by introducing a matrix A which changes 
the trigonal basis to the local spin-orbit basis in the following 
way, \a) = djj|0) = J2ia ^a,iavjcr\^)- Here wc denote the lo- 
cal spin-orbit basis using an index a(a = 1, ...6). a = 1,2 in- 
dicate the spin-orbit doublet with the total angular momentum 
Jeff =1/2 and a — 3, 4, 5, 6 denote the spin-orbit quadruplet 
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FIG. 4: Tight-binding bandstructure of Ir 5d orbitals when there 
is no trigonal crystal field effect. Here we have chosen Aso=4.0, 
ta=te=0.5. Note that the upper 4 bands corresponding to the spin- 
orbit doublet ( Jeff =1/2) are well separated from the other 8 bands 
which come from the spin-orbit quadruplet (jefi=3/l). 

with Jeff — 3/2. It is straight forward to extend the approach 
described above for a three-site cluster to the full pyrochlore 
lattice system. The resulting effective lattice hopping Hamil- 
tonian is given by 

Heff = X! d.lna^na,n'a'djn'a' + h.C, (11) 

{ij) n.n' a. a' 

where 



la.l'a' 

Here i is the unit cell index and the index n refers to the four Ir 
sites within a single tetrahedral unit cell. Including the on-site 
potentials for the local spin-orbit eigenstates, the expression 
for the effective tight-binding Hamiltonian is finally given by 
Eq.©. 

III. EVOLUTION OF THE ELECTRONIC STRUCTURE 
UNDER TRIGONAL CRYSTAL FIELD 

In this section, we discuss the evolution of the electronic 
structure of the Ir 5d system on the pyrochlore lattice in the 
presence of the trigonal crystal field. We start with the case 
where the Ir ions are under the perfect cubic crystal field with 
large spin-orbit coupling. In Fig. ID we plot the electron band 
dispersion along the high symmetry directions in the Bril- 
louin zone for the perfect cubic crystal field on the Ir ions 
with A5o=4.0, ta=te=0.5. Since we have 4 sites within a unit 
cell and each site supports 3 t2g orbitals, there are 24 bands 
within the first Brioullin zone including spin degrees of free- 
dom. Due to the time reversal and inversion symmetries, each 
band is doubly degenerate. In Fig. |4] the upper 4 bands are de- 
rived from the spin-orbit doublets with total angular momen- 
tum jeff=l/2. On the other hand, the lower 8 bands come from 



the quadruplets with jeff=3/2. Note that these two groups of 
bands are well separated by a large energy gap with an energy 
scale given by Xso- 

Since each Ir atom contributes 5 electrons, 20 bands among 
the 24 bands are filled. Namely, we have a band insulator with 
the half-filled jeff=l/2 bands. Therefore in the forthcoming 
discussion, we neglect the fully occupied jeff=3/2 bands and 
focus on the properties of the upper 4 bands (or 8 bands count- 
ing the double degeneracy of each band) possessing jeff=l/2 
character. The energy dispersions of the jeff = 1/2 states are 
shown in Fig. [Sja). The fully occupied lower two bands are 
separated from the upper two bands by a finite gap between 
them. 

To understand the topological properties of the insulating 
phase, we compute the Z2 topological invariants {i/; i'ii'2i^3) 
from the parity eigenvalues ^m(ri) at the time reversal invari- 
ant momenta, following Fu and Kane.— Here £,m{^i) indicates 
the inversion parity of the mth occupied jeff — 1/2 band at 
the time-reversal invariant momentum Fi. Using the recipro- 
cal lattice vectors {i=l, 2, 3), the eight time reversal invari- 
ant momenta can be written as ri=rn„2n3 =("-iGi + ri2G2 + 
n3G3)/2 with rii^2,3 — 0,1. The strong Z2 topological in- 
variant v is given by 

2 

(-1)'= n n ^™(r„,„,„3), (13) 

where the parity eigenvalues at the eight time reversal invari- 
ant momenta are multiplied at the same time. On the other 
hand, each of the three weak Z2 topological invariants Vi 
(i=l,2,3) is determined by the parity eigenvalues at the four 
time reversal invariant momenta lying on a plane, which is 
given by 

2 

i-^r'= n n ^™(rn,n,n3). (14) 

ni — l,nj^i— 0,1 m—1 

Because of the time reversal symmetry, each band is doubly 
degenerate at the time reversal invariant momentum and every 
Kramers doublet share the same inversion parity. Since the 
Z2 topological invariants count the parity of one state for each 
Kramers pair,'' we consider the product of the inversion par- 
ities corresponding to the two occupied jeff = 1/2 bands in 
Eq.(fT3]l and (fT4l i. Notice that, since the product of the inver- 
sion parities of the occupied jeff = 3/2 bands is H-1 in every 
time-reversal-invariant momentum, we can neglect the contri- 
butions from the jeff = 3/2 bands. These analyses lead to a 
strong topological insulator with the Z2 invariants (1;000) as 
found earlier by Pesin and Balents.^" It is interesting to note 
that a strong topological insulator with the same Z2 invariant 
(1;000) was also found in a simple one -band model on the 
pyrochlore lattice42i 

Now we describe the effect of the trigonal crystal field on 
the electronic structure of the jeff=l/2 bands. As mentioned 
above, the trigonal crystal field effect can be described by 
changing the relative magnitude of ta and t^- In general, 
the pyrochlore oxides have oxygen x parameters ranging from 
0.309 to 0.355."*" Since 2:^=0.3 125 for the perfect cubic crystal 



( a)ta/te=1.0 (b)ta/te=1.7 (c)ta/te=2.4 



(a)te/ta=1.0 (b)te/ta=2.2 (c)te/ta=2.5 
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FIG. 5: Dispersions of jcff=l/2 bands under trigonal crystal field 
witli ta > te. Here we have fixed \so ~ 4.0, te = 0.5 while 
changing ta- There is an accidental band touching between the two 
bands in the middle at the L point when ta/t^ = 1.7. 
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TABLE I: Inversion parities of the jcff=l/2 bands at time reversal 
invariant momenta for \so ~ 4.0, ta = — 0.5. Here Ei < 
E2 < Es < E4. At the momentum Xi(i=l, 2, 3), the upper (lower) 
two bands are degenerate with opposite inversion parities satisfying 
ee = -Kt'e' = -1). 



FIG. 6: Dispersions of jcff=l/2 bands under trigonal crystal field 
with ta < te. Here we have fixed Xso ~ 4.0, ta = 0.5 but vary 
tg. As te/ta increases the energy gap between the two bands in the 
middle reduces. When t^/ta ~ 2.3 band inversion occurs at the T 
point and the system becomes metallic. 
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FIG. 7: Phase diagram as a function of the strength of the lo- 
cal trigonal crystal field represented by te/ta- Note that a metallic 
phase occurs when te/ta becomes larger than the critical value of 
{te/ta)c > 1. The red cross in the middle of two strong topological 
insulators (STI) indicates the point where accidental band touching 
occurs. 



field, we have to consider both trigonal compression (x > Xc) 
and elongation (x < Xc) cases. According to the naive crys- 
tal field splitting picture, the on-site energies Ea and of 
aig and e'g states are determined by the magnitude of oxygen 
X parameters. However, in real materials, the relative magni- 
tude between Ea and E^, are strongly affected by hybridiza- 
tion with high energy eg orbital states, which are allowed un- 
der the trigonal crystal field. Therefore iiTespective of the 
magnitude of the oxygen x parameter, we have to investigate 
both ta/te > 1 and t^/ta > 1 cases on equal footing. 

We first consider the case of ta/te > 1. In Fig.|5] we plot 
the evolution of the band structure as we increase ta/t^- For 
ta/te ~ 1 .7, an accidental band touching occurs at the L point 
in the Brillouin zone. Since the two bands touching at the L 
point have opposite inversion parities as shown in Table U the 
band crossing induces the exchange of the parities between 
the two band touching at the L point. However, since we 
have eight different momentum points within the first Bril- 
louin zone, which are symmetry equivalent to the L point, the 
product of the inversion parities for all occupied bands is in- 
variant. According to Fu and Kane,^ the product of the inver- 
sion parities of occupied bands determines the Z2 topological 
invariant ly, characterizing the topological properties of insu- 
lators. Therefore the accidental band touching does not in- 
duce the change of the topological properties of the insulating 



states. 

Now we consider the opposite limit of t^/ta > 1. Fig.|6] 
shows the evolution of the band structure as we increase te/ta- 
Notice that the band gap at the F point reduces progressively 
as te/ta increases. In particular, when te/ta={te/ta)c='2--3, a 
band inversion occurs at the F point and the system becomes 
metallic. The double degeneracy at the F point is protected 
by the lattice point group symmetry. Therefore the metallic 
phase is stable as long as symmetry breaking fields reducing 
the lattice symmetry are not introduced. It is interesting that 
a similar metallic state is predicted by a recent first principle 
calculation on Y2lr207."^'* Although we have used a simplified 
tight-binding approach, the overall band structure and degen- 
eracies at the F point for the metallic phase are consistent with 
the prediction of the LDA calculation. 

The effect of the trigonal crystal fields on the electronic 
structure is summarized in the phase diagram shown in Fig.|7] 
The strong topological insulator (STI) is stable against the 
trigonal crystal field effect for te/ta smaller than the critical 
value, {te/ta)c- On the other hand, if te/ta is larger than the 
critical value, a metallic phase occurs. The metallic phase is 
stable as long as the point group symmetry of the lattice is 
protected. 
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IV. LATTICE DISTORTION AND METAL-INSULATOR 
TRANSITION 

In this section we study the fate of the metallic phase pre- 
dicted above against external perturbations. Some pyrochlore 
iridates A2lr207 with A= Nd, Sm, and Eu, show metal insula- 
tor transitions as the temperature decreases. ^^'^^ However, the 
nature of the insulating ground state is under controversy. Ac- 
cording to a recent Raman scattering measurement on these 
iridium compounds,"'^ metal-insulator transitions accompany 
structural distortions for Sm2lr207 and Eu2lr207. In addi- 
tion, a recent theoretical study on a toy model for the topo- 
logical insulators on the pyrochlore lattice,'*^ has shown that 
the lattice distortions along the [111] direction may be impor- 
tant for the realization of the topological insulator Here we 
investigate the effect of lattice distortions on the fate of the 
metallic state. In particular, to make predictions relevant to 
future high pressure experiments, we focus on uniform lattice 
modulations via q = phonon mode softening, which lead to 
macroscopic volume change of the system. 

Since we have 4 Ir ions within a tetrahedral unit cell, dis- 
placements of 4 Ir atoms from their equilibrium positions, 
Sri = {dxi,Syi,6zi), give rise to 12 independent degrees 
of freedom. Among these 12 different modes, we neglect 
6 modes describing global translations and rotations of the 
unit cell, because these do not lead to distortions of the tetra- 
hedron. The remaining 6 modes are classified as Ai sin- 
glet Q^, E doublet = (Qf , Qf ), and Ta triplet = 
{QJ, Q2, Q's), in terms of irreducible representations for the 
Td point group of the tetrahedroni^^ The expressions for 
normal coordinates Q corresponding to each irreducible rep- 
resentation are given by 







Sr23 -i 


- Sr24 -i 


- Sri3 4 


' 6r34 + 6ri2) 


Qf = 




f Sr23 


+ Sr 24 


+ Sri3 


- 25r34 - 2<5ri2) 


= 




5r2z + 


5r24 4 


Sri3) 




Qf = 




5r23) 








Qf = 




5ri3) 








Qf = 




5ri2) 






(15) 



In the above (5r,j indicates the change of the distance be- 
tween the ith and jth Ir atoms. Since the mode describes 
a uniform elongation or contraction of all bond lengths, it does 
not change the symmetry of the unit cell. So we neglect the 
Qa mode in the following discussion. On the other hand, the 
two components of the doublet = (Qf jQf) describe 
tetragonal and orthorhombic distortions, respectively. Finally, 
each component of the triplet mode, Qf depicts elongation 
and contraction of a pair of orthogonal bonds which are lying 
on two parallel planes of a cube. (See Fig. [Ub)-) Equal am- 
plitude superposition of three components of the triplet mode 




FIG. 8: Tetragonal lattice distortion induced by a q = phonon 
mode with Eg symmetry. Red solid and blue broken lines describe 
two inequivalent bond lengths under the lattice distortion. 



leads to a trigonal distortion of a tetrahedron along [1 1 1] or its 
equivalent directions. 

Using the information about the phonon modes of a sin- 
gle tetrahedron, we consider various q = phonon modes 
on the pyrochlore lattice system. The pyrochlore lattice con- 
sists of two inequivalent tetrahedra sharing a corner and these 
two types of tetrahedra are interchanged via inversion sym- 
metry with respect to a corner Therefore even when we are 
restricted to q = phonon modes, there are lattice distor- 
tions with even and odd symmetries with respect to an in- 
version center However, in this work, we focus on phonon 
modes with even inversion parities. This is because only the 
phonon modes with even parities lead to macroscopic distor- 
tion of the lattice. Therefore these modes can be softened 
via coupling to applied external pressure, which can be per- 
formed in future high pressure experiments. In particular 
we focus on tetragonal and orthorhombic distortions induced 
by Eg phonon modes and trigonal distortions driven by T2g 
phonon modes. Structures of distorted lattices and modula- 
tions in bond lengths are described in Fig. [8] and |9] 

Modulation of interatomic distances between neighboring 
Ir atoms results in the renormalization of the nearest neigh- 
bor hopping amplitudes. To understand the effect of electron- 
phonon coupling on the ground state properties of the system, 
we consider a Hamiltonian given by. 



(16) 

where ri„,7i' indicates modulation of the hopping amplitude 
between neighboring nth and ri' th Ir atoms. Specifically, we 

assume that Tna,n'a'i{Vn,n'}) =Tna,n'a'{{Vn,n' = 0})(1 - 

Vn.n')- This is equivalent to scaling of ta and te to ta{l — 
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E(k) 




FIG. 11: Band structure change induced by a trigonal lattice distor- 
tion, which comes from a q = T2g phonon mode softening. A band 
gap opens at the T point. The resultant insulating phase is a strong 
topological insulator. 



FIG. 9: Trigonal lattice distortion induced by a q = phonon mode 
with T2g symmetry. Red solid and blue broken lines describe two 
inequivalent bond lengths under the lattice distortion. 




FIG. 10: Ground state energy per unit cell for the pyrochlore iridate 
system coupled with a trigonal phonon mode along [111] direction. 
Here Aso=4.0, ta=Q-5, and te./ta=2.5. We plot the change of the 
ground state energy induced by the trigonal lattice distortion as a 
function of the hopping amplitude modulation rjT- 



rjn.n') and te(l — Vn,n') between neighboring sites n and n' . 
(See Eg. ([Toll and (fT2]i.) Since ta oc oc ctp^^^, this approx- 
imation captures the change in the overlap integral between 
neighboring d and p orbitals caused by electron-phonon cou- 
pling. The elastic constant, Kn^n', corresponding to the mod- 
ulation 77„_„', is simply taken to be Kn.n' = l-lta. 

In Fig. [To] we plot the ground state energy of the coupled 
electron-phonon system (Eq.(fT6]l) as we increase the magni- 
tude of hopping amplitude modulation rjT, which corresponds 
to a trigonal lattice distortion. Here 777- > means that the 
hopping amplitudes along the bonds connected with the site 1 
(the broken lines in Fig. |9]l are reduced by rjT while the hop- 
ping amplitudes along all other bonds are increased by the 
same amount. In Fig. (TT] we plot the change of the elec- 



Lattice distortion 


Ground state 


Tetragonal {rjEi > 0) 
Tetragonal {rjEi < 0) 


Topological Insulator 
Metal 


Orthorhombic (77^2 > 0) 
Orthorhombic (rjE2 < 0) 


Metal 
Metal 


Trigonal (rjr > 0) 
Trigonal (rjr < 0) 


Topological insulator 
Metal 



TABLE II: Lattice distortions and the resulting ground states. 



tron band dispersion induced by a trigonal lattice distortion. 
The trigonal lattice distortion results in opening a full gap at 
the Fermi energy leading to an insulating phase. The com- 
petition between the electronic energy gain from gap opening 
and elastic energy cost compromises at the equilibrium bond 
distance. Straightforward calculation of Z2 topological invari- 
ants shows that the resultant insulating ground state is a strong 
topological insulator with Z2 invariants (i'; //i 7/21^3 )=(1;000). 
Notice that the Z2 invariants of the new insulating phase are 
the same as those of the original topological insulating phase, 
which exists when there is no local trigonal crystal field split- 
ting effect. 

It is interesting that two insulating phases separated by a 
metallic phase in between, share the same topological prop- 
erties. To understand the reason for the identical topological 
properties of two insulating phases, we have checked the in- 
version parities of all bands at the T point displayed in Table|I] 
According to Fu and Kane,^ for a system with time-reversal 
and inversion symmetries, the Z2 topological invariants are 
given by the product of inversion parities of all occupied bands 
at time reversal invariant momenta. As one can see in Fig.|6] 
the metallic phase induced by local trigonal crystal field effect 
has a band crossing only at the F point. Since every band has 
even inversion parity at the F point, the band crossing does 
not change the parities of occupied bands. Therefore two in- 
sulating phases with an intervening metallic phase share the 
same Z2 topological indices in this case. 

In addition to trigonal lattice distortions, we have also in- 
vestigated the effect of tetragonal and orthorhombic distor- 
tions, which result from softening q — Q Eg phonon modes. 
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It turns out that a tetragonal distortion with rjEi > leads to 
an insulating ground state. Fig.|8]describes a bond modulation 
pattern for a tetragonal distortion with ijei > 0. Here elec- 
tron hopping amplitudes are increased for dotted bonds while 
the hopping amplitudes for solid bonds are reduced. Relative 
magnitude of bond length modulations is consistent with the 
Qf mode, the first component of the doublet E phonon mode 
in Eq.lfTsTl. However, the magnitude of the band gap generated 
by a tetragonal distortion with j]ei > is much smaller than 
that from a trigonal lattice distortion when the magnitudes of 
hopping amplitude modulations are the same in the two cases. 
In addition, an orthorhombic distortion does not open a full 
band gap. Therefore we propose that application of external 
pressure along the [1 1 1] or its equivalent directions is the most 
efficient way of generating strong topological insulators. 



V. CONCLUSION 



alent directions may lead to a strong topological insulator by 
taking the advantage of the electron-lattice coupling or per- 
haps even to a topological Mott insulator in stronger correla- 
tion regime. 

In the current work, we did not study the electron-electron 
interaction effect. It is possible that sufficiently strong elec- 
tron interaction would turn the metallic state induced by the 
trigonal crystal field effect to a magnetically ordered insulator 
or a more subtle form of Mott insulator The competition be- 
tween magnetically ordered Mott insulators, topological band 
insulators, and topological Mott insulators in the presence of 
electron interaction would be an excellent topic of future stud- 
ies. The understanding of the delicate interplay between the 
electron interaction, spin-orbit coupling, local crystal field ef- 
fect, and electron-lattice coupling would be essential for the 
determination of the ultimate ground states in these systems. 
All of these possibilities await for further experimental verifi- 
cations and findings. 



We investigate the possible existence of the strong topo- 
logical insulators in pyrochlore iridates A2lr207. The effec- 
tive Hamiltonian for the pyrochlore lattice of Ir ions is de- 
rived by taking into account the strong spin-orbit coupling 
and trigonal crystal field effect. It turns out that the strong 
topological insulator found for the ideal cubic environment of 
Ir ions may turn into a metallic state under trigonal distor- 
tion of oxygen octahedra. It has been known that, in cubic 
pyrochlore oxides, A2B2O7, where both A and B site ions re- 
side on two distinct interpenetrating pyrochlore networks, the 
trigonal crystal field splitting exists inherently '^'^ with an en- 
ergy scale comparable to the spin-orbit coupling.'*'* This may 
also be consistent with the recent LDA calculation of ¥2^207 
where it was found that the non-magnetic ground state would 
be a metal. '*'*'^*' Given that various pyrochlore iridates with dif- 
ferent A-site ions possess substantial but different amount of 
trigonal distortion, the presence of trigonal crystal field effect 
may be one of the important factors that determine the na- 
ture of low temperature ground states and finite temperature 
metal-insulator transitions. 

On the other hand, we found that the electron-lattice cou- 
pling also plays an important role. It is shown that certain q=0 
normal modes lead to the re-emergence of a strong topolog- 
ical insulator when these modes are softened and the system 
undergoes a structural deformation. Recent experimental ob- 
servation of the pressure-induced metal-insulator transition in 
Bai-^Rj^IrOs (R=Gd, Eu) suggests that electron-lattice cou- 
pling strongly affects the ground state properties in iridium 
oxide compounds. Due to the sensitive response of the elec- 
tronic structure near the Fermi level against the variation of 
Ir-O-Ir bond angles, the application of moderate hydrostatic 
pressure around 12 kbar destabilizes the metallic ground state 
leading to metal-insulator transition.'*'^ While the identifica- 
tion of the true ground state in the pyrochlore iridates would 
require better understanding of the electron correlation effect, 
it is an intriguing possibility that the strong topological insula- 
tors may arise via a finite temperature metal-insulator transi- 
tion with structural changes. Our study also suggests that the 
application of the external pressure along [111] or its equiv- 
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Appendix A: Explicit matrix expressions for the effective 
hopping Hamiltonian 

In this appendix, we present the expressions for various 
matrices used to derive the effective hopping Hamiltonian in 
Sec. Ill Bl We follow the same convention taken by Pesin and 
Balents.^° The matrix i?^*'' rotating the global cubic axis to the 
local octahedral axis for the ith Ir atom is given by 
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_ 2 


^1 






2^ 
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3 



2 2 

3 3 




2 2 
a 3 



The matrix Z)''' corresponding to a spinor representation of 
is given by 
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Finally, we show the matrix A which changes the trigonal 
basis to local spin-orbit eigenstates as follows. 



where = 3 and [nin2"-3] indicates the local C3 symme- 
try axis of the corresponding distorted oxygen octahedron. 



A = 






— 712 

6 



"3 














7llLU 

3 





\ 














—niuj 
3 







6 













-"2 






— ins 















Ve 





-ins 
V6 








/ 



^ X. -L. Qi and S. -C. Zhang, Phys. Today 63, 33 (2010). 
^ M. Z. Hasan and C. L. Kane, arXiv: 1002.3895 (unpublished). 
^ C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005). 
" C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005). 
^ B. A. Bernevig, T. L. Hughes, and S. -C. Zhang, Science 314, 
1757 (2006). 

^ B. A. Bernevig and S. -C. Zhang, Phys. Rev. Lett. 96, 106802 
(2006). 

^ L. Fu, C. L. Kane, E. J. Mele, Phys. Rev. Lett. 98, 106803 (2007). 
^ J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) (2007). 
' L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

L. Fu and C. L. Kane, Phys. Rev. B 74, 195312 (2006). 
" R. Roy, Phys. Rev. B 79, 195321 (2009). 

R. Roy, Phys. Rev. B 79, 195322 (2009). 
" M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, L. W. 

Molenkamp, X. -L. Qi, and S. -C. Zhang Science 318, 766 (2007). 

D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. 

Z. Hasan, Nature (London) 452, 970 (2008). 

Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. 

Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nat. Phys. 5, 398 

(2009). 

H. Zhang, C. -X. Liu, X. -L. Qi, X. Dai, Z. Fang, S. -C. Zhang, 
Nat. Phys. 5, 438 (2009). 
" D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F Meier, J. Oster- 
walder, L. Patthey, J. G. Checkelsky, N. P. Ong, A. V. Fedorov, H. 
Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. Hasan, 
Nature 460, 1101 (2009). 

Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. Mo, X. L. 

Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, S. C. Zhang, I. R. 

Fisher, Z. Hussain, and Z.-X. Shen, Science 325, 178 (2009). 
" A. Shitade, H. Katsura, J. Kunes, X. -L. Qi, S. -C. Zhang, and N. 

Nagaosa, Phys. Rev. Lett. 102, 125403 (2009). 
^° D. A. Pesin and L. Balents, Nat. Phys. 6, 376 (2010). 
^' B. J. Kim, H. Jin, S. J. Moon, J.-Y. Kim, B.-G. Park, C. S. Leem, 

Jaejun Yu, T. W. Noh, C. Kim, S.-J. Oh, J.-H. Park, V. Durairaj, 

G. Cao, and E. Rotenberg, Phys. Rev. Lett. 101, 076402 (2008). 
B. J. Kim, H. Ohsumi, T. Komesu, S. Sakai, T. Morita, H. Takagi, 
T. Arima, Science 323, 1329 (2009). 

S. J. Moon, H. Jin, W. S. Choi, J. S. Lee, S. S. A. Seo, J. Yu, G. 
Cao, T. W. Noh, and Y. S. Lee, Phys. Rev. B 80, 1951 10 (2009). 

H. Jin, H. Jeong, T. Ozaki, and J. Yu, Phys. Rev. B 80, 075112 
(2009). 

Y. Okamoto, M. Nohara, H. Aruga-Katori, and H. Takagi, Phys. 
Rev. Lett. 99, 137207 (2007). 



J. M. Hopkinson, S. V. Isakov, H. -Y. Kee, and Y. B. Kim, Phys. 
Rev. Lett. 99, 037201 (2007). 

M. J. Lawler, H. -Y. Kee, Y. B. Kim, and A. Vishwanath, Phys. 
Rev. Lett. 100, 227201 (2008). 

M. J. Lawler, A. Paramekanti, Y. B. Kim, and L. Balents, Phys. 
Rev. Lett. 101, 197202 (2008). 

Y. Zhou, P A. Lee, T. -K. Ng, and F. -C. Zhang, Phys. Rev. Lett. 

101, 197201 (2008). 

D. Podolsky and Y. B. Kim, |arXiv:0909.4546l (unpubUshed). 
M. R. Norman and T. Micklitz, Phys. Rev. B 81, 024428 (2010). 

G. Chen and L. Balents, Phys. Rev. B 78, 094403 (2008). 

M. W. Young, S. -S. Lee, and C. Kalhn, Phys. Rev. B 78, 125316 
(2008). 

H. Jin, H. Kim, H. Jeong, C. H. Kim, and J. Yu, arXiv:0907.0743l 
(unpublished). 

J. S. Gardner, M. J. P. Gingras, and J. E. Greedan, Rev. Mod. Phys. 
82, 53 (2010). 

K. Matsuhira, M. Wakeshima, R. Nakanishi, T. Yamada, A. Naka- 
mura, W. Kawano, S. Takagi, and Y. Hinatsu, J. Phys. Soc. Jpn. 
76, 043706 (2007). 

D. Yanagishima and Y. Maeno, J. Phys. Soc. Jpn. 70, 2880 (2001). 
T. Hasegawa, N. Ogita, K. Matsuhira, S. Takagi, M. Wakeshima, 
Y. Hinatsu, and M. Udagawa, J. Phys.: Conference. Series 200, 
012054 (2010). 

N. Taira , M. Wakeshima and Y. Hinatsu, J. Phys.: Condens. Mat- 
ter 13, 5527 (2001). 

M. A. Subramanian, G. Aravamudan, and G. V. Subba Rao, Prog. 
Sohd State Chem. 15, 55 (1983). 

M. Maczka, M. L. Sanjuan, A. F. Fuentes, L. Macalik, J. Hanuza, 
K. Matsuhira, and Z. Hiroi, Phys. Rev. B 79, 214437 (2009). 
H. -M. Guo and M. Franz, Phys. Rev. Lett. 103, 206805 (2009). 
S. Landron and M.-B. Lepetit, Phys. Rev. B 77, 125106 (2008). 
Heungsik Kim et al., to be submitted. 

M. Tinkham, Group Theory and Quantum Mechancies, (McGraw- 
Hill, New York, 1964). 

O. Tchernyshyov, R. Moessner, and S. L. Sondhi, Phys. Rev. B 

66, 064403 (2002). 

O. Tchernyshyov and G.-W. Chem, |arXiv:0907.i693l (unpub- 
lished). 

X. Wan, J. Dong, S. Y. Savrasov,| ^iv:1003.3414l (unpubUshed). 
O. B. Korneta, S. Chikara, S. Parkin, L. E. DeLong, P 
Schlottmann, and G. Cao, Phys. Rev. B 81, 045101 (2010). 



